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Abstract 

We have developed an event generator for direct-photon production in hadron collisions, in¬ 
cluding associated two-jet production in the framework of the GR@PPA event generator. The 
event generator consistently combines 7 -|- 2-jet production processes with the lowest-order 7 + 
jet and photon-radiation (fragmentation) processes from QCD 2-jet production using a subtrac¬ 
tion method. The generated events can be fed to general-purpose event generators to facilitate 
the addition of hadronization and decay simulations. Using the obtained event information, we 
can simulate photon isolation and hadron-jet reconstruction at the particle (hadron) level. The 
simulation reasonably reproduces measurement data obtained at the LHC concerning not only the 
inclusive photon spectrum, but also the correlation between the photon and jet. The simulation 
implies that the contribution of the 7 -I- 2-jet is very large, especially in low photon-pT (< 50 
GeV) regions. Discrepancies observed at low pr, although marginal, may indicate the necessity 
for the consideration of further higher-order processes. Unambiguous particle-level definition of 
the photon-isolation condition for the signal events is desired to be given explicitly in future 
measurements. 


1 Introduction 

High-energy isolated photon production in hadron collisions is considered to be suitable for probing 
the dynamics of partons (quarks and gluons) inside hadrons [2 [5], as such photons are expected 
to be produced via parton-level quantum-electrodynamic (QED) interactions, referred to as direct- 
photon production. In addition, an understanding of this process is important for hadron-collision 
experiments because various new phenomena can be explored by observing events involving energetic 
isolated photon(s). In fact, photon detection played a crucial role in the discovery of the Higgs 
boson [21S]. 

The lowest-order processes of direct-photon production are the 7 -I- jet production processes illus¬ 
trated in Fig. [ 1 ] in which the produced high-energy quark or gluon is likely to resolve to a hadron jet. 
However, these processes cannot sufficiently explain reported measurement results. A next-to-leading- 
order (NLO) calculation by JETPHOX [5] has shown that the next-order 7 -|- 2-jet contribution is 
very large; the correction to the lowest order amounts to more than 100% under typical measurement 
conditions. Typical diagrams of the 7 -I- 2 -jet production processes are shown in Fig. | 2 l Together 
with the simple correction diagrams, (a) and (b), new processes, (c) and (d), appear at the same 
coupling order. These new processes may make large contributions, as a result of the high gluon den¬ 
sity in hadrons for the process (c) and the variety of quark combinations for (d). Owing to the large 
contribution of the 7 -I- 2 -jet, NLO predictions still have relatively large energy-scale dependences. 

The CERN LHC has already provided high-statistics direct-photon data from proton-proton col¬ 
lisions to the ATLAS [21 ElISIISIUS and CMS [TTJIIIIIIS] experiments. Although the statistics is very 
high, the measurement precision is still limited. This is because the measurements suffer from the 
effects of a large hadron-decay background. As high-energy hadrons are produced in a jet-like struc¬ 
ture, this background can be dramatically reduced by imposing an isolation condition on the photon 
candidates. The isolation condition is usually defined using the sum of the transverse energies of par- 
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Figure 1: Lowest-order diagrams of direct-photon production in hadron collisions. 



Figure 2: Typical Feynman diagrams of 7 -I- 2-jet production. The diagrams (a) and (b) are simple 
gluon-radiation corrections to the diagrams in Fig. [TJ while (c) and (d) are new processes that first 
emerge at this order. 


tides inside a cone around the photon (i?™"®). The cone size (Ai?) is usually defined by the quadratic 
sum of the separations in the azimuthal angle ((/)) and pseudorapidity ( 77 ) as AR^ = -|- Arj'^, and 

E"™"® is required to be smaller than a certain value as 

( 1 ) 

Although the reduction obtained using this approach is very effective, relatively large systematic 
uncertainties remain because the remaining background events that mimic the direct-photon produc¬ 
tion are very unusual and less known. Together with the remaining energy-scale dependence in the 
NLO calculations, these uncertainties currently limit the precision of the comparison with theoretical 
predictions to the 10% level. Despite imposed large uncertainties, however, agreement between pre¬ 
diction and experimental results is marginal in the measurements reported to date, although obvious 
discrepancies have not been found. 

The measurement uncertainties will be reduced as our understanding of the background is im¬ 
proved, and theoretical predictions will become more accurate through the inclusion of further higher 
orders. However, in anticipation of comparisons with accuracy better than 10%, we have some concerns 
regarding the signal definition which is relevant to both the background and efficiency estimations. 
As we can see in the diagrams shown in Fig. [2j the 7 -I- 2-jet processes have final-state QED collinear 
divergences. The isolation conditions imposed to reduce the background also strongly reject 7 -I- 2-jet 
events in which the photon is produced almost in parallel to a quark. Those events to be rejected 
are not interesting in terms of the direct-photon measurement. In addition, inclusion of the isolation 
condition can reduce the uncertainties concerning non-perturbative effects which should regularize the 
collinear divergence. Hence, it is preferable to include the isolation condition in the definition of signal 
events for which measurement results and theoretical predictions are presented. 

In order to conduct rigorous comparisons between measurements and predictions, we must define 
the isolation condition unambiguously. The detection efficiency must be estimated for events satisfying 
the isolation condition, and those events that do not satisfy the isolation condition must be treated as 
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background, even if they are produced via one of the processes shown in Fig. [2l Usually, theoretical 
calculations apply isolation conditions at the parton level. However, in principle, it is impossible to give 
unambiguous definitions at the parton level since partons are not physical objects. The appearance 
of partons depends on the perturbation order and arbitrary cutoffs. Unambiguous definitions can be 
given only by using particle(hadron)-level information, as in the hadron-jet definition |14] . Such a 
definition has already been adopted in recent measurements Ellin], in which the particle species to 
be included for the evaluation of U™"® are also explicitly specified. 

If the definition is given at the particle level, we must therefore convert parton-level theoretical 
predictions into predictions at the particle level. Although it may be possible to evaluate relevant 
corrections using separate simulations, the most straightforward approach is to use parton-level Monte 
Carlo (MC) event generators that can be connected to simulations down to the particle level. In such 
simulations, it is preferable that the event generation is based on matrix elements (MEs) consistently 
including multi-jet production processes, in order to achieve sufficient precision. 

We have been developing event generators that facilitate the matching between processes having 
different jet-multiplicities, based on a subtraction method ESI ES] in the framework of the GR@PPA 
event generator EZHU]. The generated events can be fed to general-purpose event generators, such as 
PYTHIA ES] and HERWIG ED] , in order to simulate hadronization and decays to obtain particle-level 
event information. We first developed event generators that combine weak-boson production processes 
associated with 0 -jet and 1 -jet production miiia, and demonstrated that the developed event gen¬ 
erator reproduces measured Z-boson px spectra very precisely ElEllSSj- In these event generators, 
collinear divergent components of quantum-chromodynamics (QGD) are subtracted from 1-jet MEs in 
order to prevent double-counting with the jet production in parton-shower (PS) simulations applied to 
0-jet events. The matching method has been extended to diphoton ( 77 ) production [DD], in which the 
subtraction is applied to both QED and QGD collinear divergences. Recently, we have succeeded in 
combining 2 -jet production processes with 0 - and 1 -jet production where the subtraction method 
is modified to account for the soft-gluon divergence together with the collinear divergence. 

In this article, we construct an event generator that consistently combines the simulation of 7 -I- 
2 -jet production with the lowest-order 7 -|- 1 -jet using the techniques that we have developed to date. 
QGD 2-jet production having hard-photon radiation from PS is also incorporated in order to restore 
the subtracted QED collinear components. The generated events, to which further simulations down 
to the particle level are added, are compared with measurement results at the LHG. Although the event 
generation is based on tree-level MEs only, we can expect the simulations to have precision comparable 
to NLO calculations, as the correction to the lowest order is expected to be dominated by 7 -I- 2-jet 
real radiation contributions. Subsequently, we present other distributions that can be derived from the 
simulation results, in order to discuss the characteristic properties of direct-photon production. The 
understanding of this process obtained through such studies will be helpful in reducing the systematic 
uncertainties in future measurements. 

Note that a similar event generator has already been developed in the SHERPA framework pS] , 
and its simulation results reproduce the event kinematics of diphoton ( 77 ) |29j and direct-photon [Dj 
production with good precision. In the SHERPA event generator, processes having different jet- 
multiplicities are merged using the GKKW method m, in which multi-jet events generated according 
to jet-including MEs are suppressed by reinterpreting them in the context of a PS. The suppression 
is applied without separating non-divergent components from divergent components that PS is based 
on. Such an ME-based simulation is applied down to a small merging scale, and application of PS 
is limited to further soft regions. We are taking a different approach to achieving the same goal by 
using the subtraction method. In our approach, the subtraction and PS simulations are applied up to 
the typical energy scale of events for which MEs are evaluated. Thus, no suppression is required and 
non-divergent components in MEs are preserved as they are in the entire phase space. 

The rest of this article is organized as follows. The overall simulation strategy is discussed in 
Sec. m and the self-consistency of the developed event generator is examined in Sec. [Dj We compare 
the simulation results with LHG measurements in Sec.|4| Using the simulated events, the characteristic 
properties of direct-photon production are discussed in Sec. [51 and the discussions are concluded in 
Sec. ini 
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2 Simulation strategy 


Programs for calculating the MEs of the 7 + 2-jet production processes were newly generated using 
the GRACE system [JT] [32] and installed in the GR@PPA event generator framework. The MEs are 
based on Eeynman diagrams including all possible combinations of quarks (up to the b quark) and 
gluons having photon radiation from one of the quark lines. Thus, the QCD part of the process is 
identical to the full set of QCD 2-jet production processes from which gg —>• gg is excluded. 

The 7 -I- 2-jet MEs have various collinear divergences. We subtract these divergences using the 
limited leading-logarithmic (LLL) subtraction method [H] to render the cross sections finite. Together 
with the QCD divergences, final-state QED divergences are subtracted using a method that has been 
developed previously m- We do not need to consider initial-state QED divergences since an energetic 
photon is required to be detected at large angles. In the subtraction of QCD divergences, we also need 
to account for the soft-gluon divergence in the 7 -|- 2-jet production. This divergence is subtracted 
simultaneously with the collinear divergences, using the combined subtraction method developed in 
our previous study m- 

The subtraction is applied in order to prevent double-counting of the jet production in MEs and 
that from PSs attached to lower jet-multiplicity events. Thus, the subtracted components are restored 
by combining lower multiplicity processes to which PS simulations are applied. The subtracted QCD 
components are restored by combining 7 -|- 1-jet production processes. In accordance with the modi¬ 
fication to the subtraction designed to account for the soft-gluon divergence, the 7 -|- 1 -jet events are 
reweighted after PS simulations are applied EZj. 

The subtracted QED components are restored by combining QCD 2-jet production processes to 
which a QCD/QED-mixed PS is applied [33]. The photon radiation in low-Q^ regions, which PS 
simulations do not accommodate, is simulated by using a fragmentation function (PE) [33] as in the 
diphoton-production simulation m- In order to improve the generation efficiency, the production of 
one energetic photon is enforced in this PS/FF simulation. The enforcing procedure determines an 
additional event weight for the simulated events. 

The event generation is performed using the LabCut framework of GR@PPA m- In this frame¬ 
work, PS simulations are applied before differential cross section values are passed to the integration 
and event generation utility, BASES/SPRING [331133] . employed in GR@PPA. Thus, the cross section 
values can be weighted using event weights that are determined after application of the PS simula¬ 
tions. In the generation of 7 -|- 1-jet and QCD 2-jet events, the evaluated event weights are accounted 
for in the integration and event generation using this framework. 

The LabCut framework was originally introduced for the efficient generation of events for which 
strict kinematical constraints are required. As PS simulations alter the momenta of generated particles, 
the constraints in parton-level hard interactions must be relaxed significantly in order to avoid biases 
to the final sample. Direct-photon production is such a process. For instance, when we generate the 
events used in the next section, in which a pr cut of 100 GeV/c is applied to the photon, we must 
reduce the threshold to 40 GeV/c in the hard-interaction generation. We would suffer from extremely 
poor efficiency in the event selection if we adopted an ordinary generation scheme. For this reason, 
we also adopt the LabCut framework for the generation of 7 -I- 2-jet events. 

Since PS simulations are limited by a certain value (Q^ < Mpg), the subtraction must also be 
limited by the same in order to achieve matching between them. However, setting the limit for 
the LLL subtraction (/iLLu) equal to pps is not trivial, since these values are determined in different 
processes. In order to solve this problem, GR@PPA includes a mechanism that refers to energy scales 
defined for different processes. In the LLL subtraction, /ipg of the non-radiative event is referred to 
using event information that is reconstructed by removing the assumed radiation, in order to assign 
the same value to /illl- Thus, in principle, any energy-scale definition can be adopted to achieve 
matching. 

The above energy scale may differ for the initial-state radiation (ISR) and the final-state radiation 
(FSR). Hence, we denote these boundaries as pisn and pfsb., respectively. In the event generation, 
we must also specify the other energy scales necessary for QCD calculations, i.e., the renormalization 
scale (pr) and the factorization scale (pf)- As we require a basic identity between the initial-state 
PS and the QCD evolution in parton distribution functions (PDFs), pisn must be equal to /Xi? in our 
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matching method. In principle, fin may differ from However, since Og running is assumed in PS 
and PDF, it would be better to choose an identical value in order to achieve an equal jet-production 
strength at the boundary between PS and ME. For these reasons, we use an identical value for fip, 
fip, and /iisR as the default setting. Although there is no strict reason, we also adopt the same value 
for ^fsr in the present study, for simplicity. Thus, the energy scale is unique in the present study, 
j.e., fj, = fiR = Up = fiisR = fJ-FSR- 

We must then decide how to define /r. As we have discussed above, any definition is acceptable for 
the matching. However, a uniform definition that is applicable to all considered processes is preferable. 
Direct-photon production is considered to be dominated by t-channel interactions. Hence, it would 
be natural to use the pp value of the photon for fi. However, the photon is not necessarily the hardest 
particle in 7 -|- 2-jet production. Moreover, no photon appears in QCD 2-jet production. 

In the present study, we adopt a new definition that can be applied to all considered processes, 

fi = max{(5T,i}- (2) 

The quantity QT,i is defined for each final-state particle in terms of its mass and pp, such that 

QT,i = ^^+Ph^ (3) 

and the largest value is taken as the energy scale of the event, as defined in Eq. ©. This definition 
reduces to p = m for single-resonance production and to = rn? + p\ for pair production. It also 
reduces to p^ = my -\-pp for vector boson {V) + jet production processes. In the processes considered 
in the present study, p is determined by the pp value of the particle having the largest pp since all 
particles are nearly massless. 

The generated event samples of 7 -I- 2-jet (a2j), 7 + 1-jet (alj), and QCD 2-jet {qcd2j), to which 
PS simulations are fully applied down to Qo = 5 GeV, are passed to PYTHIA 6.425 [H] in order to 
add hadronization and decay simulations. Although low-Q^ (< Qq) QCD effects are also simulated 
in PYTHIA, their effects are almost negligible in terms of the observable quantities of direct-photon 
events. However, the momenta of the generated particles are more or less altered by these simulations. 
In order to avoid possible biases, the generation conditions for the events to be passed to PYTHIA are 
slightly relaxed from the final selection criteria. Por instance, accounting for a large safety margin, the 
Pp threshold of photons is reduced to 90 GeV/c for the final selection condition of pp > 100 GeV/c. 
The PYTHIA simulation is applied with its default setting, except for the settings of PARP(67) = 
1.0 and PARP (71) = 1.0 as in the previous studies El [Ml US]. In addition, the underlying-event 
simulation is deactivated (MSTP(81) = 0) since LHC experiments correct the isolation cone energy 
for this effect. 

As described earlier, the PS simulations in GR@PPA rotate and boost the original events gener¬ 
ated by hard interactions significantly. In order to check possible biases due to the hard-interaction 
generation conditions, the event information prior to application of the PS simulations is passed to 
PYTHIA, together with the PS-applied event information. The spectra relevant to the original gener¬ 
ation conditions are checked after applying the particle-level event selection, in order to confirm that 
the selected sample is not biased by the generation conditions. 


3 Matching test 


The self-consistency of the developed event generator was tested by investigating observable quantities 
in the simulated events. The events were generated under the 7-TeV LHC condition, i.e., for proton- 
proton collisions at a center-of-mass (CM) energy of 7 TeV, using the CTEQ6L1 PDF [31] included 
in the GR@PPA package. We required that the events simulated down to the particle level contain 
the photon generated by GR@PPA in the transverse-momentum (pp) and pseudorapidity ( 7 ) ranges 
of 


Pt{j) > 100 GeV/c and 17 ( 7 )| < 2.4. 


(4) 


As particle-level event information is available, we can apply a realistic isolation condition to the 
photons. We adopted the simplest definition in Eq. O- The cone Ep (A™"®) was evaluated using all 
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stable particles, excluding neutrinos and muons, produced inside the cone around the photon having 
a size of Ai? = 0.4. We defined the isolation condition as 

< 7.0 GeV. (5) 

We need hadron-jet information in the simulated events in order to test the matching of the 7 -I- 
2-jet production with other processes. The hadron jets were reconstructed using Fast Jet 3.0.3 [37], 
with the application of the anti-^T algorithm with R — 0.4. All stable particles in the pseudorapidity 
region of I 77 I < 5.0, including neutrinos, were used for the reconstruction. The reconstructed jets that 
satished the conditions, 

PT(jet) > 30 GeV/c, | 77 (jet)| < 4.4, and Ai?( 7 ,jet) > 0.5, ( 6 ) 

were taken as the detected jets to be used in the subsequent analyses, where Ai?( 7 , jet) is the separation 
in Ai?, Ai?^ = A(/)^ -I- IS.rf', between the photon and jet. We required that at least two jets were 
detected in each event. 

Since ^ defines the boundaries between the PS and ME-based simulations for divergent compo¬ 
nents, a strict test of the matching can be conducted by investigating the stability of observable 
distributions against the variation of fji. We take the scale defined by Eq. m as the standard value 
/tq. The simulations were carried out for extreme choices of /r = /ro/2 and /i = 2/ro, together with the 
standard setting of /r = /tq. The pr threshold of photons in Eq. 0 was selected in order to ensure 
that pL is always larger than the pr threshold of hadron jets in Eq. 0 . Hence, events are always 
allowed to include soft hadron jets that are predominantly produced by PS simulations. 

Because of a large contribution from the 7 -I- 2-jet, the integrated cross section has a significant 
/i dependence, yielding a 37% increase for p = po/2 and an 18% decrease for p = 2po with respect 
to the value obtained for the standard setting. Nearly two-thirds of the variation is caused by the 
change in pn and the remainder is due to pp- As these changes are not relevant to the matching to 
be tested here, we compare the simulation results in terms of the normalized differential cross section, 
{lja)dafdx, where x denotes the tested quantity. 

The matching for the final-state photon radiation can be tested by investigating the correlation 
between the photon and the second jet, i.e., the detected hadron jet having the second-largest px- 
The distribution of the Ai? separation between them is shown in Fig.[3| The qcd2j result and the sum 
of the alj and a2j results, labeled alj + a2j, are shown separately in the upper panels. The photon is 
produced by PS in the former case, while it is produced in the hard interaction in the latter. Events 
in the Ai? < 0.5 region are rejected by the requirement imposed on the jets in Eq. ®. The p = po 
results are illustrated using solid histograms, and those for p = po/2 and p = 2pQ are represented by 
dashed and dotted histograms, respectively. We can see that these two results have remarkably large 
p dependences. The sum of the results for all simulated processes is shown in the lower panel. The 
large p dependences observed in the separate results are compensated by each other to form a very 
stable distribution in the summed result, implying that almost complete matching is accomplished. 

Analogously, Ai? separation between the leading (highest-pr) and second jets is a good measure 
to test the matching in the hnal-state QCD radiation. This Ai? distribution is shown separately for 
alj and a2j -|- qcd2j in the upper panels in Fig. |4| The second jet is produced by PS in the former, 
while the hard interaction contains two jets in the latter. The sum of these results is shown in the 
lower panel. The events are not distributed in the Ai? < 0.4 region because of the cone size in the 
jet reconstruction. Again, although the dominant ISR contributions to the second jet render the FSR 
contributions less clear, we can see that the large p dependences observed in the separate processes 
are compensated to yield a stable distribution in the summed result, implying good matching in the 
FSR. 

Despite an overall compensation, a small p dependence is left in the summed result in Fig.U] near 
the lower limit of Ai?. This is due to the multiple radiation effects in PS. Our matching method takes 
care about the leading contribution of PS only. The method is tailored so as to achieve matching 
between this leading contribution and the divergent terms in radiative processes. The matching may 
be deteriorated by the existence of additional jets in PS in the range where p^ is moved. On the 
other hand, only one photon is assumed in the present study because the contribution from multiple 
photon radiation is expected to be negligible. This must be the reason why we observe almost perfect 
compensation in Fig.[3| 
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Figure 3: Distribution of the AR separation between the photon and second jet. The distributions 
of qcd2j and alj + a2j are shown separately in the upper panels, and the sum is shown in the lower 
panel. The results for /i = /io/2, /io, and 2/2o are represented by dashed, solid, and dotted histograms, 
respectively. 
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Figure 4: Distribution of the AR separation between the leading and second jets. The distributions 
for alj and a2j + qcd2j are shown separately in the upper panels, and the sum is shown in the lower 
panel. The results for n = /xo/2, /iq, and 2/io are represented by dashed, solid, and dotted histograms, 
respectively. 
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Figure 5: pr distribution of the second jet. The distributions for alj and a2j + qcd2j are shown 
separately in the upper panels, and the sum is shown in the lower panel. The results for p = /io/2, 
POi and 2/io are represented by dashed, solid, and dotted histograms, respectively. 


The matching in the initial-state QCD radiation can be tested by investigating the second-jet pr 
distribution. The alj and a2j -|- qcd2j results are shown separately in the upper panels in Fig. [Sj and 
the sum is shown in the lower panel. We again observe a strong compensation for the p dependence in 
this result, implying good matching in the initial-state QCD radiation. However, some p dependence 
remains in the summed result, although this may be difficult to see in Fig. [5] where the cross sections 
are displayed on a logarithmic scale. The reason for the remaining dependence must be identical to 
that discussed for the results shown in Fig. |4l 

In conclusion, we have observed a strong compensation between separate processes for the energy- 
scale dependence of the observable quantity distributions, which are sensitive to the matching in the 
QED final-state and QCD final- and initial-state radiations. These observations imply good matching 
between photon and jet radiation in PS and those in MEs of radiative processes. We can obtain stable 
distributions by adding the results of all processes. The remaining energy-scale dependence can be 
attributed to a multiple-radiation effect in PS. 


4 Comparison with measurements at LHC 

The performance of the developed event generator is examined in this section by comparing the 
simulation results with actual measurement data obtained at the LHC. Because a large amount of 
measurement data is currently available, we compare our results with selected measurements below. 


4.1 ATLAS measurement (1) 

We neglect to consider the first CMS m and ATLAS [HI measurements because they are based on 
limited data obtained at an early stage of these experiments. In this subsection, we compare our 
simulation results with the second ATLAS measurement [7] based on 35-pb“^ data obtained at a CM 
energy of 7 TeV and collected in 2010. They measured the photon Et (= Pt) spectra in the range of 
45 < Et < 400 GeV in four pseudorapidity (ry) regions: \r]\ < 0.6, 0.6 < |? 7 | < 1.37, 1.52 < \r]\ < 1.8, 
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Figure 6: Comparison with the ATLAS measurement data reported in [7]. The simulation results are 
represented by boxes. The vertical sizes of the boxes indicate the statistical error of the simulation. 
The measurement results are shown with plots having error bars, which indicate the total error quoted 
in Tables A.l - A.4 in the ATLAS paper, where the luminosity uncertainty of 3.4% is not included. 
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Figure 7: Ratio of the simulation results to the ATLAS measurement data reported in [7]. The 
notation is identical to that in Fig. 
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and 1.8 < I 77 I < 2.37. They applied the cone isolation condition in Eq. ([T]) with E^° = 3 GeV at the 
detector level, and estimated that this condition is equivalent to EJf° = 4 GeV with Ai? = 0.4 at the 
particle level, in which all stable particles except for neutrinos and muons are included. 

Following the measurement conditions, we simulated the events in 7-TeV proton-proton collisions, 
with the conditions Et > 45 GeV, \ri\ < 2.37, and i?™"® < 4 GeV, and binned them to derive 
the differential cross section. A modified leading-order PDF, MRST2007LO* [38], was used for the 
hard-interaction generation because this PDF is widely used in experiments for MC simulations. The 
generation conditions and preselection conditions for the events to be fed to PYTHIA were sufficiently 
relaxed in order not to bias the simulation results. Effectively, no constraint was applied to associated 
jets. 

The simulation results are compared with the measurement data in Fig. ( 6 ] and the ratio between 
them is presented in Fig. |7| Because the measurement error is dominated by systematic errors, the 
simulation was carried out separately for \ow-Et (< 100 GeV) and high-Ex (> 100 GeV) regions in 
order to save the computational cost. Although there is no guarantee of the overall normalization for 
simulations based on tree-level MEs, the simulation results are in good agreement with the measure¬ 
ment data within the quoted measurement errors, without any adjustment of the normalization. 

Despite an overall agreement, we can observe a slight enhancement of the simulation in Iow-Et 
regions. The energy scale /i is another free parameter in the simulation that we can optimize. However, 
a change of /r does not result in significant alteration of the Et dependence, although it alters the 
overall yield. Even if we change the value from /io/2 to 2^o, as in the study of the matching in Sec. (3] 
the variation in the normalized cross section, ja^da jdEx■, is confined within ± 10 % from the fj, = fj,o 
result. The use of a leading-order PDF, GTEQ6L1 [35], reduces the overall yield by approximately 
10%, but it does not significantly alter the Ex dependence. This must be simply due to the lower 
gluon density in this PDF. 

4.2 CMS measurement (1) 

The inclusive direct-photon cross section was also measured by GMS [12] based on 36-pb“^ data 
collected in 2010. They measured the differential cross section, d^a/dEx/dt], in the Ex range of 25 
< Ex < 400 GeV in four pseudorapidity ranges within \ri\ < 2.5. Although an isolation condition 
is required in the measurement, the signal definition is not explicitly described. As they apply the 
condition, A™"® < 5 GeV with Ai? = 0.4, to the simulated events in their discussions, we assume this 
condition to be the signal definition in the present study. 

We carried out a simulation by appropriately setting the Ex threshold and the 77 range, and applied 
the above isolation condition to the simulated events. The other simulation conditions are identical to 
those in Sec. 01 The A™''® was evaluated using all stable particles except for neutrinos and muons, 
although this detail is not given in the CMS paper. Again, in order to save the computational cost, 
we performed the simulation in three separate Ex ranges: 25 < Ex < 50 GeV, 50 < Ex < 100 GeV, 
and 100 <Ex< 400 GeV. 

We present the ratio between the simulation and measurement in Fig. [5] Although the simulation 
is in good agreement with the measurements at high Ex 100 GeV), a significant excess of the 
simulation is observed in lowev-Ex regions. This tendency is consistent with the observation in 
Sec. 14.11 CMS compared their results with an NLO prediction by JETPHOX [5] and observed a 
similar excess of the prediction at low Ex- However, the excess that they observed is smaller than 
the present one. As we will show later, the 2-jet contributions {a2j ± qcd2j) increase as Ex decreases. 
Hence, the contribution of further higher-order processes may also be large at low Ex- JETPHOX 
partly involves 3-jet contributions as an NLO correction to the fragmentation processes. The observed 
excess of the simulation may be related to the contribution from missing higher-order processes. 

4.3 ATLAS measurement (2) 

ATLAS measured the properties of photons, jets, and their correlations in 7 ± jet events [S]- The 
measurement was based on the same data as the measurement discussed in Sec. sn with the same 
photon selection condition. In order to study the jet properties, they reconstructed hadron jets using 
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Figure 8 : Comparison with the CMS measurement data for inclusive direct-photon cross section 
reported in m- The ratio of the simulation to the measurement is illustrated using boxes. The 
vertical size of each box represents the statistical error of the simulation. The error bars indicate the 
quadratic sum of the statistical and systematic errors of the measurement presented in Table 6 of the 
CMS paper, where the luminosity uncertainty of 4.0% is not included. 


the anti-ZcT clustering with R = 0.6, and applied the conditions, 

PT(jet) > 40 GeV/c, |r?(jet)| < 2.37, Ai 7 ( 7 ,jet) > 1.0, (7) 

to the reconstructed jets. The notation for these parameters is identical to that in Sec. [21 

In order to simulate the measurement results, we applied the same jet clustering to the simulation 
events used in Sec. 14. II employing Fast Jet 3.0.3, and also applied the kinematical condition in Eq. m 
to the reconstructed jets. All stable particles in the pseudorapidity range of \r]\ < 5.0, including 
neutrinos, were used for the clustering. The observation of at least one such jet was required, and 
the leading (highest-pr) jet was taken as the jet to be studied. The simulation results are compared 
with the measurements in Fig. |9l The presented results correspond to those in Figs. 9-14 and 17, 
and Tables 1 - 7 in the ATLAS paper [9] . The plotted error bars represent the quadratic sums of the 
statistical and systematic errors listed in the tables. 

The agreement between the simulation and measurement is marginal. The majority of the discrep¬ 
ancies that we observe in Fig. [5] are caused by the small excess of the simulation results at low At ( 7 ) 
and low pT(jet). We observe nearly complete agreement for all results, excluding the A(/)( 7 , jet) distri¬ 
bution, if we assume an additional 20% inefficiency around pT(jet) = 50 GeV/c for the jet detection, 
or 10% photon inefficiency plus 10% jet inefficiency for At ( 7 ) and pT(jet) around 50 GeV. 

The discrepancy in the A(/)( 7 , jet) distribution, i.e., the distribution of the azimuthal angle separa¬ 
tion between the photon and jet, cannot be significantly reduced by assuming additional inefficiencies. 
As described in Sec. |2l we deactivated the underlying-event simulation in PYTHIA. Angular distri¬ 
butions are sometimes very sensitive to such soft interactions. However, we observed no significant 
improvement when the underlying-event simulation was activated. Rather, this simply reduced the 
overall yield because of an additional contribution to A™"® for the photon isolation. 

A similar discrepancy in the A^( 7 ,jet) distribution was observed by ATLAS when their measure¬ 
ment result was compared with predictions from JETPHOX and HERWIG, while the SHERPA and 
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Figure 9: Comparison with the ATLAS measurements for 7 + jet events reported in [^. The simula¬ 
tion results are represented by boxes and the measurement results are plotted with error bars. The 
luminosity uncertainty of 3.4% is included in the error bars. 
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PYTHIA predictions exhibited better agreement (see Fig. 12 in i)- SHERPA includes multi-jet ME 
contributions, and PYTHIA simulates multiple radiation effects in the collinear approximation by 
increasing the PS energy scale. Hence, the discrepancy in A(^( 7 ,jet) observed in the present study 
must be due to the lack of multi-jet (three or more jet) contributions in our simulation. 

4.4 CMS measurement (2) 

CMS also measured the properties of 7 -I- jet events in terms of the triple-differential cross section, 
a/{dp^dr]jdr]jet), using 2.14-fb“^ data for 7-TeV collisions [13], in the pseudorapidity range of 
|? 7 | < 2.5 for both photons and jets, and in the pr ranges of 40 < pt < 300 GeV/c for photons and 
Pt > 30 GeV/c for jets. The hadron jets were reconstructed using the anti-fey clustering with R = 
0.5, and the separation between the photon and jet was required to be Ai?( 7 , jet) > 0.5. 

The simulation was carried out following the above conditions by using the MRST2007LO* PDF. 
Although the isolation condition for the photon signal is again not explicitly defined, we assumed the 
condition if™"®(Ai? = 0.4) < 5 GeV in our simulation because CMS required this condition in their 
simulation. In order to save the computational cost, the simulation was carried out separately for the 
Pt{i) ranges below and above 85 GeV/c. 

The ratio of the simulation to the measurement data is presented in Fig. (TUI We observe reasonable 
agreement between them; at least, no remarkable systematic discrepancy is observed. A marginal 
observed discrepancy is a systematic dehcit of the simulation in the angular region of 1.57 < 17 ( 7 ) | < 
2.1 and | 7 (jet)| < 1.5. However, a deficit is observed in the normalization only and the discrepancy is 
marginal compared to the measurement errors. We also observe a marked discrepancy in the highest- 
Pt{j) bins in the largest-] 7 ( 7 )! regions. This discrepancy was also observed by CMS in comparison 
with the JETPHOX prediction in their report. We have no comment on this observation. We must 
note that, in this result, we do not observe the excess in the simulation result at low pt that we 
observed in Sec. 

4.5 ATLAS measurement (3) 

ATLAS published inclusive Et{'^) distribution results in high-Fly regions using 4.6-fb“^ data collected 
in 2011 [TU] ■ They presented the differential cross section, da / dEx , in the Ay range of 100 GeV < Ay < 
1 TeV in two pseudorapidity ranges, \r]\ < 1.37 and 1.52 < jp] < 2.37. The definition of the isolation 
condition for the signal photons was explicitly given as A|P"®(Ai? = 0.4) < 7 GeV, where A™"® was 
evaluated by using all stable particles included in the cone, excluding neutrinos and muons. The 
simulation was carried out following these measurement conditions using the MRST2007LO* PDF. 
The events were simulated separately for Ay ranges below and above 250 GeV. 

The ratio of the simulation results to the measurement data is presented in Fig.[TTJ Gontrary to the 
previous comparison results, the simulation exhibits a systematic deficit of approximately 10 % through 
the measurement range. We observe very good agreement with the measurement if we increase the 
overall normalization of the simulation by 10%. As nothing is essentially changed in the simulation, 
this discrepancy indicates an inconsistency between the measurements. However, since the error of 
the measurements discussed previously was at the level of 10 % or larger, this inconsistency could be 
considered as a tolerance within the measurement errors. 

4.6 Summary of the comparison 

Overall, the simulation based on the developed MC event generator reproduces direct-photon mea¬ 
surement results at the LHG reasonably well, not only for inclusive photon spectra but also for the 
properties of associated jet production. Here, we emphasize that this jet is not an energetic parton, 
but is a realistic hadron jet reconstructed from simulated stable particles. 

Despite an overall consistency, some discrepancies seem to exist in low-py regions (< 50 GeV/c). 
This observation may indicate the necessity for further higher-order processes including three or more 
jets in the final state. However, this observation is not conclusive as the measurements are not fully 
consistent with each other. New measurements or reanalyses are necessary for confirmation. 
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Figure 10: Comparison with the CMS measurement data for 7 + jet events reported in [T3]. The ratio 
between the simulation and measurement concerning the triple-differential cross section is presented 
as a function of the photon px- The results are presented in the eight angular regions, four 177 ( 7 )! ^ 
two |77(jet)|, shown in the figure. The error bars include the luminosity uncertainty of 2.2%. 
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Figure 11: Comparison with the ATLAS measurement data for high-Ay ( 7 ) regions published in | 10 j . 
The ratio of the simulation to the measurement concerning the inclusive cross section is represented 
by boxes as a function of At ( 7 ) in two angular regions. The error bars indicate the quadratic sum 
of the statistical and systematic errors of the measurement, including the luminosity uncertainty of 
1 . 8 %. 


As we have mentioned in the introduction, the signal events for which the measurement results 
are presented must be unambiguously defined to permit rigorous comparisons with theoretical pre¬ 
dictions. However, the definition of the photon isolation condition is not explicitly presented in some 
measurement reports. We would like to find clear descriptions in future reports. In addition, we 
are concerned about the estimation of the photon-isolation efficiency. The efficiency is evaluated in 
the majority of the measurements using the PYTHIA simulation. As we have noted in a previous 
study [26], the photon-radiation probability is too small in the PYTHIA 6.4 PS simulation. We are 
concerned that this problem may have somewhat affected the measurement results published to date. 


5 Characteristic properties 

In this section, we discuss the characteristic properties of the direct-photon production process in 
hadron collisions using our simulation results. We first examine the impact of the photon-isolation 
requirement. We can evaluate the cross section without applying the isolation cut, as the cross section 
is made finite by summing all-order contributions employing PS and FF in collinear regions. Figure[T2] 
shows the efficiency of the isolation cut, be., the ratio between the cross sections with and without 
the cut. This simulation was carried out under the 7-TeV LHC condition using the MRST2007LO* 
PDF with our standard energy-scale setting. The direct photon was required to be produced in the 
angular region of \ri\ < 2.4. 

The efficiency was evaluated for three values of the simple cone-A t cut defined in Eq. CD with 
AA = 0.4 as a function of the photon pT- It can be seen that the efficiency is approximately 80% and 
decreases slightly as pt{i) increases. It almost saturates at high pt{i) 200 GeV). The dependence 
on the AT-cut value is moderate in the examined range. The simulation shows that approximately one 
half of the qcd2j events are rejected by the cut. On the other hand, only a few percent of the alj and 
a2j events are rejected, because these processes do not have any enhancement in photon-quark collinear 
regions. Although the rejection by the isolation cut is not very large, it is inadvisable to correct the 
measured cross section for this inefficiency, because substantial uncertainties are suspected in the 
efficiency evaluation, such as uncertainties in non-perturbative effects and PS-model dependencies. 

Figure |T3| shows the fractional contribution of the lowest-order 7 -|- 1 -jet production processes 
after the isolation cut of A“"® < 7 GeV. Here, we evaluated the ratio of the 7 -I- 1-jet cross section 
without applying the soft-gluon correction with respect to the summed cross section. The result is 
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Figure 12: Efficiency of the photon-isolation cut. The ratio between the cross sections with and 
without the cut is displayed as a function of the photon pT- The simulation results under the 7-TeV 
LHC condition are presented for three cut values of = 4, 7, and 15 GeV. 
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Figure 13: Fraction of the lowest-order 7 -|- 1-jet contribution after the isolation cut of E™''® < 7 GeV 
with our standard energy-scale setting. The ratio of the 7 -I- 1-jet cross section without the soft-gluon 
correction to the summed cross section, a\j -I- a2j + qcd2j, is displayed as a function of the photon px- 
The simulation results obtained by using two different PDFs, MRST2007LO* (solid) and CTEQ6L1 
(dashed), are presented. 
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Figure 14: Distribution of the AR separation between the photon and second jet. The results for alj 
and a2j + qcd2j are shown separately. Three results for energy scales of p, = po/2 (dashed), p = po 
(solid), and p = 2po (dotted) are presented. The summed distribution with p = po is illustrated using 
thick dots for comparison. 


presented as a function of pt{i)- The soft-gluon correction reduces the 7 -I- 1-jet cross section, but 
the reduction is not larger than 10% through the examined pt{i) range. We can see that the fraction 
is only approximately 50% at pt{i) = 30 GeV, despite the fact that a large portion of qcd2j events 
have been rejected by the isolation cut. In other words, as we have emphasized, the correction by the 
7 -I- 2-jet contributions is quite large, being nearly 100%. Hence, we suspect that further higher-order 
contributions may also be large at low pt[i)- In Fig.lT^ the results of the simulations conducted using 
two PDFs, MRST2007LO* and CTEQ6L1, are presented. The difference between the two results is 
only a few percent. 

Note that the result in Fig. [T3] is not physically meaningful. The sharing between the processes 
depends on the choice of energy scales. However, as we adopted a natural choice of energy scales, the 
above hndings must be helpful for furthering understanding of direct-photon production. In theoretical 
calculations, the choice of p = pT {'j)/2 is also frequently adopted. The lowest-order fraction is further 
reduced for this choice. 

As we have shown in Sec. [31 the final-state QED divergence produces a flat distribution of events in 
the distribution of the AR separation between the photon and second jet, AR{j, 2nd jet), at small AR 
(^ 1-5)- The existence of final-state divergences is characteristic in higher-order processes involving 
two or more jets. The contributions of the lowest-order 1-jet (alj) and higher-order 2-jet {a2j + qcd2j) 
processes to this AR distribution are shown separately in Fig. [Ml We used the simulation sample for 
the matching test in Sec. [31 and applied the same analyses to the events. The only difference is in the 
separation of the processes. 

As in Sec. [31 the results for the energy-scale choices of p = po/2 and p = 2po are shown together 
with those for the standard choice of p = po. In this result, we observe a very large energy-scale 
dependence in large-Ai? regions. However, the dependence reduces as AR decreases, and the flat 
distribution at AR <1.5 cannot be reproduced by the lowest-order alj contribution. Thus, if we can 
measure this distribution at small AR, this will be a direct quantitative measurement of higher-order 
contributions that necessarily have QED divergences in the final state. 

Finally, we examine the correlation between the photon and the leading jet. As the direct-photon 
production is a two-body 7 -I- jet production at the lowest order, this process is considered to constitute 
a good sample for jet-energy calibration; the pt of hadron jets can be calibrated by referring to 
the photon px which can be measured precisely. However, higher-order processes may deteriorate 
this capability. The distribution of the ratio between the photon px and the leading-jet px, Rx = 
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Figure 15: Distribution of the pt ratio between the photon and leading jet, Rt = PT(lst jet)/p 7 ’( 7 )- 
The simulation sample used in the matching test in Sec. [3] was used here. The photon and jets were 
selected according to the conditions in Eqs. O, ®, and requiring the detection of at least one 
jet. The results for four pt{j) ranges are presented. 





Figure 16: Identical to Fig. [131 but the number of detected jets is required to be exactly one. 
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Priest jet)/pT{7), in the events used for the matching test discussed in Sec. [3] are shown in Fig. [Tsl 
The photon and jet were selected according to the conditions in Eqs. and and at least 

one jet was required to be detected. As a result of the very loose condition for jet detection, at least 
one jet is almost always detected in photon-detected events; the jet-detection efhciency is more than 
99%. 

The Rt distribution is shown separately for four pripf) ranges in the figure. Although a concen¬ 
tration of events is observed around Rt = 1, we observe a shoulder-like structure at small Rt- This 
structure is created by the contribution of second jets produced in association with the leading jet. 
Thus, it is expected that this structure can be eliminated by suppressing the second-jet contribution. 
As an example, we required that the number of detected jets was exactly one. Approximately one half 
of the events were rejected by this requirement. The loose jet-detection condition became a tight veto 
condition in this part of the study. Figure [TCI shows the results obtained under this requirement, and it 
is apparent that the structure at small Rt has vanished. In addition, the long tail to the Rt > 1 side 
has been eliminated. This is because this tail was predominantly produced by qcd2j contributions, in 
which a detectable jet was produced in association with the photon. 

The Rt distributions in Fig.[T6]still have finite widths because of the small-angle undetectable jets 
and soft QCD interactions. In addition, especially in low pt( 7) ranges, the distribution has a longer 
tail to the Rt < 1 side. This tail is due to the particles escaping from the jet reconstruction. The 
peak position is somewhat shifted from unity, possibly as a result of this effect. 

These observations imply that direct-photon production provides a good data sample for the 
calibration of hadron-jet momenta if we can efficiently eliminate the contributions from additional 
jets. In practical applications, we must account for the shift in the peak position caused by the 
unavoidable escape of particles from jets. The shift would be larger than that in Fig. [THlin real-world 
cases, because of the additional detection resolution. 


6 Conclusion 

We have developed an event generator for direct-photon production in hadron collisions, consistently 
combining 7 -|- 2-jet production processes with the lowest-order 7 -|- I-jet production using a sub¬ 
traction method. The final-state QED collinear divergences are also subtracted from 7 -|- 2-jet MEs. 
The subtracted QED components are restored by combining QCD 2-jet production events in which 
an energetic photon is produced by a QCD/QED-mixed PS and a simulation employing an EE. The 
photon radiation is enforced in this PS/EE simulation to facilitate efficient event generation. The 
generated events can be passed to PYTHIA in order to add hadronization and decay simulations. 
We can conduct realistic simulations of photon isolation and hadron-jet clustering using the obtained 
particle-level event information. 

The self-consistency of the developed event generator has been examined by investigating ob¬ 
servable quantities in the simulated events. The obtained distributions, which are sensitive to the 
matching between 7 -|- 2-jet production and other contributions, exhibit good stability against vari¬ 
ation of energy scales, implying good matching between the included PS and ME-based simulations 
for photon and jet radiation. 

The performance of the simulation has been tested by comparing the simulation results with 
measurement data at the LHC. The comparison was carried out concerning not only the inclusive 
properties of the photon, but also the properties of the associated hadron jet. The simulation results 
show reasonable agreement with the measurement data, especially in high photon-pr regions, where 
Pt{i) ^ 100 GeV. Although substantial discrepancies have been found in low-pr regions of the photon 
and jet, ie., pt ^ 50 GeV, the observations are not conclusive because the measurement results are 
not fully consistent with each other. 

We have demonstrated some characteristic properties of direct-photon production using the simu¬ 
lated events. The efficiency of the photon-isolation cut is approximately 80% if a simple cone-Ay cut is 
applied with cut values of approximately 10 GeV. Although the cut predominantly rejects final-state 
QED radiation events, the contribution of the 2-jet processes is very large even after the isolation 
cut. The lowest-order 7 -|- 1-jet contribution is as small as 50% around pt (7) = 50 GeV with our 
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standard energy-scale setting. This fraction gradually increases as pril) increases. Hence, discrepan¬ 
cies between the simulation and measurement observed at low px may indicate the necessity for the 
incorporation of further higher-order processes, including the production of three or more jets. The 
higher-order contributions can be quantitatively measured by investigating the angular separation 
between the photon and the second jet. We have also shown that the direct-photon process provides 
a good sample for the jet-momentum calibration, if the second-jet contribution can be effectively 
eliminated. 

As we have noted in the introduction, it is important to give an unambiguous definition of the 
photon isolation for the signal events in measurements, in order to facilitate rigorous comparison with 
theoretical predictions. The isolation must be defined using the event information at the particle 
(hadron) level. We could not find such a definition for some of the measurements that we used for the 
comparison conducted in the present study. Future measurements are desired to present the definition 
explicitly. 
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